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Standard Nearest Neighbor Discretizations of Klein-Gordon Models 
Cannot Preserve Both Energy and Linear Momentum 
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We consider nonlinear Klein-Gordon wave equations and illustrate that standard discretizations 
thereof (involving nearest neighbors) may preserve either standardly defined linear momentum or 
total energy but not both. This has a variety of intriguing implications for the "non-potential" dis- 
cretizations that preserve only the linear momentum, such as the self-accelerating or self-decelerating 
motion of coherent structures such as discrete kinks in these nonlinear lattices. 
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Introduction. In the last two decades, the interplay 
of nonlinearity and spatial discreteness has been increas- 
ingly recognized as vital for the understanding of a va- 
riety of physical systems 0. Such contexts range from 
calcium waves in living cells 0] to the propagation of ac- 
tion potentials through the cardiac tissue Q and from 
chains of chemical reactions £j to applications in super- 
conductivity and Josephson junctions [j|, nonlinear op- 
tics and fiber/waveguide arrays |6|, complex electronic 
materials 0, Bose-Einstein condensates la] or the local 
denaturation of the DNA double strand |J 

On the other hand, spatially discrete systems (of cou- 
pled nonlinear ordinary differential equations) are also 
relevant as discretizations and computational implemen- 
tations of the corresponding continuum field theories that 
are applicable to a variety of contexts such as statistical 
mechanics 0, solid state physics [ill, fluid mechanics 
[l2| and particle physics 01 ( see a l so references therein). 
Nonlinear Klein-Gordon type equations are a prototyp- 
ical example among such wave models and their vari- 
ants span a variety of applications including Josephson 
junctions in superconductivity, cosmic domain walls in 
cosmology, elementary particles in particle physics and 
denaturation bubbles in the DNA, among others. 

In this communication, we examine some of the key 
properties that ensue when discretizing nonlinear Klein- 
Gordon (KG) equations, using nearest neighbor approxi- 
mations (which are the most standard ones implemented 
in the literature; see e.g., 0). In particular, we focus 
on the physically relevant invariances of the continuum 
equation (more specifically, the conservation of the linear 
momentum and of the total energy of the system) and il- 
lustrate the surprising result that if we demand that the 
energy be conserved, then the momentum cannot be con- 
served, while if we demand that the momentum be con- 
served then the energy cannot be conserved (resulting in 
a so-called non-potential model |lj|). Our presentation 
will be structured as follows. First, we will provide the 
general mathematical setting of KG equations and study 
their discretizations that conserve linear momentum and 
energy, comparing and constructing the properties of the 
two. Then, we are going to give an application of our 



considerations to the physically relevant tp 4 model, i.e., 
the ubiquitous double well potential. Finally, we will 
summarize our conclusions and discuss future directions. 

Setup and Analytical Results. We consider the La- 
grangian of the Klein- Gordon field, 
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4i - - v(4>) 



dx . 



and the corresponding equation of motion, 

4>u = (f>xx - v'{4>) . 



(i) 



(2) 



Assuming that the background potential V((f>) can be 
expanded in Taylor series we write 



(3) 



For brevity, when possible, we will use the notations 

4> n -i =1, <p n = m, cj> n+ i =r. (4) 

We start with a general proof of our main statement, 
namely that discretizations that preserve linear momen- 
tum and energy are mutually exclusive for nearest neigh- 
bor discretizations. As was shown in [l^, the standard 
discretization of Eq. J3J) that preserves the discrete ana- 
log of the linear momentum, defined in a standard way, 



M = 4> n {<t>n+l 



»n-l) 



is one of the form: 



m = C(l + r- 2m) - 



F{r,m) - F(m,l) 



(5) 



(6) 



where C = l/h 2 , where h is the lattice spacing, and 
the derivative of F is equal to V in the continuum limit 
(C -> oo). Then, 
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'M+l 







= Yjfi^n+lAn) - H{4> n , <j^_{j\ = 0, 



2 



where H(r, m) — C(r 2 + m 2 — 2mr) — F(r, m), and the 
terms 4>n{4>n+i — 4>n-i) cancel out as the telescopic sum. 

However, if the model is potential, for nearest neigh- 
bor discretizations the nonlinear term will be of the form 
V(r, m) such that the Lagrangian can be written as: 
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K+l ~ 4>n) 2 ~ V((j> n+1 ,(j> n ) 



> (8) 



where the first term gives the kinetic energy, K,, and the 
two other terms give negative potential energy, —V, so 
that the total energy is E = KL + V . However, then a 
model that would enforce both energy and momentum 
conservation would have to satisfy: 



F{r, m) - F(m, I) d 
r — I dm 



V(r, m) + V(m, I) 



(9) 



After multiplying with r — I, this, in turn, implies that 
the cross terms involving all 3 of r, m and I should be 
presentable in the form 

dV(m,l) ,dV(r,m) , , , x 

r ^—^--l ^ — - = P{r,m) - P{m,l). (10) 



dm 



dm 



This is satisfied only if V (x, y) is a (symmetric) quadratic 
function in its arguments. However, this is incompatible 
with the nonlinear nature of the model. Hence, it is not 
possible to satisfy both conservation laws at once. 

Let us now derive the general discrete Klein-Gordon 
model of the form of Eq. © conserving momentum. For 
the polynomial background forces Eq. (j3J), the nonlinear 
term of Eq. JSJ can be presented as the sum of s-order 
terms 

B(l, m ,r) = F ^-J^ =±B 8 (l, m ,r), (11) 



s=0 



with 



where 



i—0 j—i 



E E & «> = ,t - 

i—0 j=i 



(12) 



(13) 



In the continuum limit one has / — > m and r — * m and 
thus, Eq. Ijllll together with Eq. (|13f) ensure the desired 
limit, V'(4>). Furthermore, Eq. (|12|) takes into account 
all possible combinations of powers of Z, m, and r. Coeffi- 
cients bij. s make a triangular matrix of size (s+1) x (s+1). 
Let us find the coefficients 6y iS to satisfy Eq. (|llfl . Wc 
write 



('• /:/''- EE'''-''' 7 

*=0 j=i 
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i=0 j=i 



Terms containing both I and r should be canceled out 
because they do not fit the representation of Eq. lllfl . 
This can be achieved by setting 6y, s — 6(j_|_i)(j_|_i) )S , i.e., 
coefficients in each diagonal of the triangular matrix must 
be equal. The simplified expression reads: 



{r-l)B s = J2b 



, 5 s T 771 



(15) 



To symmetrize the result, we add and subtract boo, s m s+1 
(r - l)B s = b oAr s+1 + ^ s+1 ) - 6 o, s (m s+1 + l s+1 ) 

S S 

+ £ b ois _ i+1)!S r i m s - i+1 - Yl b i,smH s - i+l {lQ) 

i=l i=l 

where we shifted the summation index by 1 in the first 
sum and also used the equality of the diagonal coeffi- 
cients. The desired representation is obtained for arbi- 
trary 6 00iS and arbitrary b 0hs = 6o(«-i+i),s for i > 0. 
Summing up, (i) the coefficients bij. s within each diago- 
nal are equal, (ii) the coefficients on the main diagonal 
can be chosen arbitrarily, and (iii) the terms on zth super- 
diagonal (i > 0) must have the same coefficients as the 
terms on (s — i + l)th diagonal (and these can also be 
chosen arbitrarily). 

For B s the number of super-diagonals is s so that the 
number of free coefficients is 1+ \s /2] , where [x] is lowest 
integer greater than or equal to x. We must also take 
into account the relation between coefficients Eq. i|13|) 
and the number of free coefficients becomes \s/2]. 

For example, the coefficients of B3 are 



&01.3 
bn 



J00,3 
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'00,3 



&02,3 &01,3 
&01,3 &02,3 
&00 : 3 &01,3 
^00,3 

4&01,3 + 26q2, 3 = C3. 



(17) 



Since the model Eq. © is translationally invariant, the 
static kink is free of the Peierls-Nabarro potential (PNp) 
|15|. i.e., the periodic potential that nonlinear waves have 
to overcome to move by one lattice site [l6| (see also refer- 
ences therein). This is an important qualitative difference 
with respect to the conventional discretization when, in 
Eq. QJ, V(4>) is substituted with V(4> n ) and thus, in the 
equation of motion Eq. J2|, V((f>) becomes V'(<j> n ). 

Another class of Klein-Gordon models which support 
energy conservation and sustain static kinks but which 
are free of PNp has been derived by Speight and collab- 
orators I n such models the background potential 
term of Eq. (JTJ, P = J V((f>)dx, should be discretized as 



G{<t> n+1 )-G{4> n )\ w . th [GW] 2 = F(0) 



Un+1 



(18) 

Numerical Results/Model Comparison. We now exam- 
ine various models proposed as discretizations of the con- 
tinuum field theory in the context of perhaps one of the 
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most famous such examples, namely the double-well 4> 4 
model [HIHEi (see also the review p^. 

The discrete Klein-Gordon model conserving momen- 
tum is given by Eq. © with the nonlinear term Eq. 
where the coefficients &jj, s are as described in the previous 
section. The continuum <fr model has the background 
potential V{4>) = (1 - </> 2 ) 2 /4, hence V'{4>) = ~4> + </> 3 so 
that in Eq. J3J all a s — except for o\ = —03 = — 1. 
The momentum-preserving discretization then reads: 

to = (C + a) (I + r — 2m) + m 
-f3(l 2 + lr + r 2 ) + 0m(l + r + m 
-7 (l 3 + r 3 + l 2 r + lr 2 ) - 5m (l 2 + to 2 + r 2 + Ir) 

- 4 7 - 4<5)m 2 (l + r), (19) 

where a = -&oo,i, ^ = ^00.2, 7 = ^00,3, <5 = 601,3 are free 
parameters and we did not include the terms with s > 3. 

The model of Eq. 1|19[) will be compared to the model 
obtained from Eq. (|18|l in </> 4 case [lj], namely 

m = + (I + r - 2m) + m 

[2m 3 + (m + /) 3 + (m + r) 3 ] , (20) 
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and also to the "standard" </> 4 discretization, i.e., 

m = C(l + r -2m) + m-m 3 . (21) 

If in Eq. a = /3 = 7 = <5 = 0, then the models 

of Eq. l|19f) and Eq. (|20|l have the same linear vibra- 
tion spectrum (i.e., dispersion relation u — u>(k)) for 
the vacuum solution 4> n = ±1, namely uj 2 = 2 + (AC — 
2) sin 2 («;/2). This can be compared to the spectrum of 
the vacuum ofEq. lj?T|) . ui 2 = 2 + 4Csin 2 (W2). 

We analyze the kink internal modes (i.e., internal de- 
grees of freedom ^^) for these three models. First, we 
determine the kink-like hcteroclinic solution by means of 
relaxational dynamics. Then, the linearized equations 
are used in a lattice of N — 200 sites to obtain N eigen- 
frequencies and the corresponding eigenmodes. We are 
particularly interested in the eigenfrequencies which lie 
outside the linear vibration band of vacuum solution and 
thus are associated with the kink internal modes. It is 
worthwhile to notice that the eigenproblem for models 
conserving energy, Eq. I|20ll and Eq. (|21|l . has a symmet- 
ric Hessian matrix while the non-self-adjoint problem for 
the momentum conserving model Eq. (|19fl results in a 
non-symmetric matrix. 

The top panels of Fig. nP resen t the boundaries of the 
linear vibration spectrum of the vacuum (solid lines) and 
the kink internal modes (dots) as the functions of lattice 
spacing h for (a) the classical 4 model of Eq. H21|) , (b) 
the PNp-free model of Eq. (|20|l conserving energy, and 
(c) the PNp-free model of Eq. 119fl conserving momen- 
tum. In PNp-free models kinks possess a zero frequency, 
Goldstone translational mode similarly to the continuum 
cj) 4 kink. Hence, the static kink can be centered anywhere 
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FIG. 1: Upper panels: boundaries of the linear spectrum of 
the vacuum (solid lines) and kink internal mode frequencies 
(dots) as functions of the lattice spacing h = 1/VC. Lower 
panels: time evolution of kink velocity for different initial 
velocities and h = 0.7. The results are shown for (a) classical 
4> A model, Eq. (b) PNp-free model conserving energy, 

Eq. (J20J, and (c) PNp-free model conserving momentum, Eq. 
with a = /3 = 7 = (5 = 0. 
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FIG. 2: Trajectories of particles (a) in the model of Eq. I|19|l 
with h = 0.7 when the kink moves with a steady velocity v* 
(see Fig. ^c), bottom panel) and (b) for the continuum rj) 4, 
kink. 



on the lattice. The results presented in Fig. ^are for the 
kink situated exactly between two lattice sites. This po- 
sition is the stable position for the classical c/> 4 discrete 
kink . Since all 3 discrete models share the same con- 
tinuum (</> 4 ) limit, their spectra are very close for small 
h(< 0.5). We found that the model Eq. (O may have 
kink internal modes lying above the spectrum of vacuum, 
e.g., for a = 1/2, (3 = 0, 7 = 1/4, and (5 = 0. Such 
modes are short-wavelength ones, with large amplitudes 
(energies) and they do not radiate because of the absence 
of coupling to the linear phonon spectrum. 

Perhaps more interesting are the implications of such 
discretizations on the mobility of these localized coher- 
ent structures. In the PNp-free models, Eq. (|19|l and 
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Eq. (HI, the kink was launched using a perturbation 
along the Goldstone mode to provide the initial kick. In 
the classical model Eq. lf!HTl for this purpose we employed 
the imaginary frequency (real eigenvalue) unstable eigen- 
mode for a kink initialized at the unstable position (a 
"site-centered" kink). In all cases the amplitude of the 
mode is related to the initial velocity of the kink. In 
the bottom panels of Fig. we present the time evo- 
lution of the kink velocity for different initial velocities 
and h = 0.7 for the three discretizations. The results 
suggest that the mobility of the kink in the classical </> 4 
model presented in (a) is much smaller than in the PNp- 
free models, (b) and (c). Furthermore, a very interesting 
effect of kink self- acceleration can be observed in panel 
(c). Here there exists a selected kink velocity v* ss 0.637 
and kinks launched with v > v* , in a very short time ad- 
just their velocities to v* . More surprisingly, the veloc- 
ity adjustment is observed even for kinks launched with 
v < v* . In the steady-state regime, when the kink moves 
with v = v* , it excites (in its tail) the short-wave oscilla- 
tory mode even though in front of the kink the vacuum 
is not perturbed. 

These results generate the question of where the energy 
for the self-acceleration and vacuum excitation comes 
from. In Fig. 2(a) we show the trajectories of four 
neighboring particles when a kink moving with v = v* 
(see Fig. dc), bottom panel) moves through. For com- 
parison, in (b) the trajectories for the classical 4> A kink, 
4> n (t) = tnnh[p(nh — vt)], where p = [2 — 2v 2 }~ 1 / 2 , are 
shown. In both cases the trajectories are identical and 
shifted with respect to each other by t = h/v, but in (b) 
they are the odd functions with respect the point (j> n = 
while in (a) they are not. The work done by the back- 
ground forces, Eq. Qllfl. to move the nth particle from 
one energy well to another is W n = — J_ rhB(l, m, r)dt. 
For the 4> A model Eq. O with (3 = 7 = S = 0, 
the nonlinear part of B(l,m,r) reduces to B(l,m,r) — 
(l/2)m 2 (l + r). It is straightforward to demonstrate 
that W n = for the classical <fi 4 kink. However, if a 
term breaking odd symmetry, e.g., e cosh^ 1 [9(nh — vt)], 
is added to the kink, the work becomes nonzero, W n = 
^e(e 2 + 1) [cosh(ph) — l] 3 / sinh 4 (/?/i), where we set for 



simplicity 9 — p. Numerically we found that W n can be 
positive or negative depending on p, 9 and the kink ve- 
locity, v. This simple analysis qualitatively explains the 
kink self-acceleration or deceleration and the vacuum ex- 
citation. The energy for this comes from the breaking of 
the odd symmetry of particle trajectories, which is pos- 
sible in the case of path-dependent background forces. 
It is, thus, very interesting to highlight the distinctions 
between the "regular" discrete models, the PNp-free, en- 
ergy conserving discrete models, and the PNp-free, mo- 
mentum conserving discrete models. The first ones lead 
to rapid dissipation of the wave's kinetic energy due to 
the PN barrier. The second render the dissipation far 
slower in time. Finally the third may even sustain self- 
accelerating waves and locking to a particular speed due 
to the non-potential nature of the relevant model. 

Conclusions and Future Challenges. The statement 
that Klein-Gordon discrete model cannot conserve en- 
ergy and momentum simultaneously was proved for the 
case of standard nearest-neighbor discretizations. This 
raises the issue of how to gauge the "adequacy" of a dis- 
cretization scheme with respect to the continuum model 
dynamics. In view of this question, a number of charac- 
teristic similarities and differences between energy- and 
momentum-conserving discrete models were highlighted. 
The momentum conserving Klein-Gordon system with 
non-potential background forces discussed here differs 
from other path-dependent systems, e.g., having friction 
and/or AC drive, in the sense that the viscosity and ex- 
ternal forces are not explicitly introduced. This makes 
the dynamics of the system peculiar, for example, as it 
was demonstrated, the existence, the intensity, and the 
sign of energy exchange with the surroundings depends 
on the symmetry and other characteristics of the motion. 
Further investigation of the intriguing dynamic proper- 
ties of such non-potential models is important, given the 
relevance of path-dependent forces in various applications 
such as, e.g., aerodynamic and hydrodynamic forces, the 
forces induced in automatic control systems and others. 
Such studies are in progress and will be reported in future 
publications. 



[1] S. Flach and C. R. Willis, Phys. Rep. 295, 181 (1998); 
D. Hennig and G. Tsironis, Phys. Rep. 307, 333 (1999); 
P. G. Kevrekidis, K. O. Rasmussen, and A. R. Bishop, 
Int. J. Mod. Phys. B 15, 2833 (2001); Yu. S. Kivshar 
and G. P. Agrawal, Optical Solitons: From Fibers to Pho- 
tonic Crystals, Academic Press (San Diego, 2003); D. N. 
Christodoulides, F. Lederer, and Y. Silberberg, Nature 
424, 817 (2003). 

[2] S. P. Dawson, J. Keizer, and J. E. Pearson, Proc. Natl. 
Acad. Sci. U.S.A. 96, 6060 (1999). 

[3] J. P. Keener, J. Theor. Biol. 148, 49 (1991). 

[4] J. P. Laplante and T. Erneux, J. Phys. Chem. 96, 4931 
(1992). 



[5] A. V. Ustinov, T. Doderer, I. V. Vernik, N. F. Pedersen, 
R. P. Huebener, and V. A. Oboznov, Physica D 68, 41 
(1993). 

[6] D. N. Christodoulides and R. I. Joseph, Opt. Lett. 13, 
794 (1988). 

[7] B. L. Swanson et al, Phys. Rev. Lett. 82, 3288 (1999). 
[8] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86, 2353 
(2001). 

[9] M. Peyrard and A. R. Bishop, Phys. Rev. Lett. 62, 2755 
(1989). 

[10] G. Parisi, Statistical Field Theory, Addison- Wesley (New 
York, 1988). 

[11] R. K. Dodd, J. C. Eilbeck, J. D. Gibbon, and H. C. Mor- 



■5 



ris, Sohtons and Nonlinear Wave Equations, Academic 
Press (London, 1982). 
[12] E. Infeld and G. Rowlands, Nonlinear Waves, Solitons 
and Chaos, Cambridge University Press (Cambridge, 
1990). 

[13] P. Anninos, S. Oliveira, and R. A. Matzner, Phys. Rev. 

D 44, 1147 (1991). 
[14] V. V. Bolotin, Nonconservative problems of the theory of 

elastic stability, Pergamon Press, (Oxford, 1963). 
[15] P. G. Kevrekidis, Physica D 183, 68 (2003). 
[16] Yu. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 



763 (1989). 

[17] J. M. Speight and R. S. Ward, Nonlinearity 7, 475 
(1994); J. M. Speight, Nonlinearity 10, 1615 (1997). J. M. 
Speight, Nonlinearity 12, 1373 (1999). 

[18] T. I. Belova and A. E. Kudryavtsev, Phys. Uspekhi 40, 
359 (1997). 

[19] O. M. Braun, Yu. S. Kivshar, and M. Peyrard, Phys. 
Rev. E 56, 6050 (1997); Yu. S. Kivshar et al, Phys. Rev. 
Lett. 80, 5032 (1998); P. G. Kevrekidis and C. K. R. T. 
Jones, Phys. Rev. E 61, 3114 (2000). 



